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ABSTRACT 

We present the first measurement of the planet frequency beyond the "snow 
hne" , for the planet-to-star mass-ratio interval —4.5 < log < —2, corresponding 
to the range of ice giants to gas giants. We find 



- (0.36 ±0.15) dex" 



d log q d log s 

at mean mass ratio g = 5 x 10^"* with no discernible deviation from a fiat (Opik's 
Law) distribution in log projected separation s. The determination is based on 
a sample of 6 planets detected from intensive follow-up observations of high- 
magnification [A > 200) microlensing events during 2005-2008. The sampled 
host stars have a typical mass Mhost ~ O.5M0, and detection is sensitive to 
planets over a range of planet-star projected separations (s~^i?E, Smax-RE), where 
Re ~ 3.5AU(Mhost/M0)i/2 jg Einstein radius and s^ax ~ (g/10-^-3)2/3. 
This corresponds to deprojected separations roughly 3 times the "snow line". 
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Despite the frenetic nature of these observations, we show that they have the 
properties of a "controlled experiment" , which is what permits measurement of 
absolute planet frequency. High-magnification events are rare, but the survey- 
plus- followup high-magnification channel is very efficient: half of all high-mag 
events were successfully monitored and half of these yielded planet detections. 
The planet frequency derived from microlensing is a factor 7 larger than the 
one derived from Doppler studies at factor ~ 25 smaller star- planet separations 
[i.e., periods 2 - 2000 days]. However, this difference is basically consistent with 
the gradient derived from Doppler studies (when extrapolated well beyond the 
separations from which it is measured). This suggests a universal separation 
distribution across 2 dex in planet-star separation, 2 dex in mass ratio, and 
0.3 dex in host mass. Finally, if all planetary systems were "analogs" of the 
Solar System, our sample would have yielded 18.2 planets (11.4 "Jupiters", 6.4 
"Saturns", 0.3 "Uranuses", 0.2 "Neptunes") including 6.1 systems with 2 or more 
planet detections. This compares to 6 planets including one two-planet system in 
the actual sample, implying a first estimate of 1/6 for the frequency of solar-like 
systems. 

Subject headings: gravitational lensing - planetary systems 

1. Introduction 

To date, 10 microlensing-planet discoveries have been published, which permits, at least 
in principle, a measurement of planet parameter distribution functions. Of course, the size 
of this sample is small, both absolutely and relative to the dozens of planets that have been 
discovered from transit surveys and the hundreds from Doppler (radial velocity, hereafter 
RV) surveys. However, microlensing probes a substantially different part of parameter space 
from these other methods. The majority of RV planets, and the overwhelming majority of 
transiting planets are believed to have reached their present locations, generally well within 
the "snow line" , by migrating by a factor > 10 inward from their birthplace. By contrast, 
microlensing planets are generally found beyond the snow line, where gas giants (analogs of 
Jupiter and Saturn) and ice giants (analogs of Uranus and Neptune) are thought to form. 
Thus it would be of substantial interest to compare the properties of microlensing planets, 
which have not suffered major migrations, to other planets that have. 

In the present incarnations of microlensing experiments, planet detection occurs through 
several different routes, with various degrees of "human intervention" . At one extreme, plan- 
ets can be detected simply on the basis of survey data, without even the knowledge that a 
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planet was present until after the event is over. This is the plan for so-called second gener- 
ation microlensing planet surveys, which will per mit rigorous, particle -phys ics-like objective 



analys is of a "controlled experiment". In fact, iTsapras et al.l ( l2003l ) and ISnodgrass et al. 



(120041 ) have already carried out such analyses of first-generation survey-only data. How- 
ever, to date only one se cure microlensing pla net has been discovered from survey-only data, 
MOA-2007-BLG-192Lb JSennett et al.lboosh . 



Rather, most microlensing planet detections have taken place through a complex inter- 
play of survey and followup observations. Only a year after the first microlensing events were 
discovered, the MACHO collaboration issued an lAU circular urging follow-up observations 
of a mic rolensing event, an d the OGLE c ollaboration initiated its Early Warning System 
(EWS) fludalski et allll994j : ludalskilboosh . which regularly issued "alerts" of ongoing mi- 
crolensing events, usually well before peak, as soon as they were reliably det ected. EWS, 



toget her with a similar program soon initiated by the MACHO collaboration ( lAlcock et al. 



19961 ) with their much large r camera, enabled th e formation in 1995 of t he first microlensiri g 
follow-up teams, PLANET dAlbrow et al.lll998h and MACHO/GMAN jAlcock et al.lll997h . 
and soon thereafter, MPS ( Rhie et al.l 1999 ). The EROS team issu ed only a few a lerts, but 



one of these led to the first mass measurement of an unseen obi ect flAn et al 



first spatially resolved high-resolution spectrum of another star fICastro et al 



2002), and the 
200ll ). OGLE- 



II inaugurated wide-field observations as well in 1998. In order to cover large areas of the sky 
(even with wide-field cameras), the survey teams would generally obtain only ~ 1 point per 
night per field. Since typical planetary deviations from "normal" (point lens) microlensing 
light curves las t only of order a day or less, this was not generally adequate to detect a 
planet. Hence, iGould Sz Loebl (119921 ) already advocated the formation of follow-up teams 
that would choose several favorable events to monitor more frequently, using telescopes on 
several continents to permit 24-hour coverage. In 2000, MACHO ceased operations, but a 
new survey group, MOA, had already begun survey observations. 

Once this synergy between survey and follow-up teams was established, it evolved 
quickly on both sides. Both types of teams developed the capacity for "internal alerts", 
whereby real-time photometry was quickly analyzed to find hints of an anomaly. If these 
hints were regarded as sufficiently interesting, they would trigger additional observations by 
the team. The first such alert was by MACHO/GMAN in 1996 and several followed the next 
year from PLANET. In 2003, OGLE developed the Early EWS, which automatically alerted 
the observer to possible anomalies, who then would make additional observations and, if 
these were confirming, publicize them to the community. OGLE also pioneered making their 
data publicly available, which greatly facilitated follow-up work. These developments gener- 
ally evolved into a system of mutual alerts, open transfer of data between teams, and active 
ongoing email discussion of developing events. 
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The first secure microlensing planet, OGLE-2003-BLG-235Lb ( iBond et al.l l2004j ) was 
discovered by means of such an internal alert. The MOA team noticed deviations in their 
data, and initiated intensive follow-up of this event. In retrospect, one can see that even 
without this internal alert, the combined OGLE and MOA data would have been sufficient 
to show that the star had a companion. However, the normal survey data would not have 
permitted one to tell whether this companion was a planet or a low-mass star (or possibly a 
brown dwarf). 

Detections of this type are at the opposite extreme from the pure-survey detection, 
which can be modeled as a controlled experiment. Without understanding the efficiency at 
which the observers issue their real time alerts, one cannot measure the absolute rate of 
planets from observations that ar e triggered by th e presence of the planet itself. Another 



planet, OGLE-2007-BLG-368Lb (ISumi et al.l 1201011 falls par tlv into this category. In this 
case, the alert was triggered by ARTEMIS JPominik et aPboOsh . whose ultimate goal is 
to communicate such alerts directly to the followup telescopes without human interference 
(and so bring such alerts into the fold of "rigorous controlled experiments"), although at this 
stage ARTEMIS alerts were still vetted by humans. A further anomaly was soon detected 
by the MOA observer. But even without followup observations triggered by these alerts, 
this companion probably could have been constrained to be planetary, although with larger 
errors, so this detection could still be integrated into the "controlled experiment" framework 
even without trying to model the alert process. 

However, if follow-up observations are carried out without being significantly infiuenced 
by the possible presence of a planet, then these also can be treate d as a controlled e xperi - 
men t, and absolute ra tes can be derived using either the method of iGaudi fc SackettI (120001 ) 
or of iRhie et al.l (120001 ). Such an analysis was car r ied out for the firs t 43 events monitored by 
the PLANET Collaboration (lAlbrow et al.ll200ll : iGaudi et al.ll2002l ). In particular, since no 
planets were detected (and so only upper limits obtained), there was no possibility of recog- 
nition of a possible planet playing any role in the observation strategy. But, as emphasized 
above, once planets started to be discovered, the situation became more nuanced. 

One path toward obtaining followup observations that are triggered on potentially plan- 
etary anomalies, without compromising the "controlled experiment" ideal, is to establish 
robotic triggering programs that communicate directly to robotic telescopes. This has been 
the goal of the RoboNet Collaboration since its inception in 2004 jTsapras et al.ll2009l). Al- 
gorit hms for detecting anomalies in rob otically acquired re al-time data ( iDominik et al.ll2007l . 
20081 ) and directing robotic telescopes ( iHorne et al.ll2009l ) have been devised and are being 
further developed, and will expand their scope as the Las Cumbres Observatory Global 
Telescope Network (LCOGT) itself expands. But to date, the most effective RoboNET 
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observations have been "hand triggered" (e.g., for OGLE-2007-BLG-349Lb). 



In spite of this nuanced situation, ISumi et al.l ( 120101 ) were nevertheless able to extract 



relative frequencies of planets from the ensemble of all 10 microlensing planets. They argued 
that each of these planet detections was characterized by a sensitivity g{q) oc g™, where 
q is the planet/star mass ratio. While m varies from event to event, they argued that 
m = 0.6 ±0.1 was an appropriate average value. They then fit the ensemble of detections to 
a power-law mass distribution f{q)d\ogq oc q'^dlogq and found n = —0.68 ± 0.20. But they 
did not attempt to extract absolute frequencies from their sample. 

Here, we analyze an important subclass of microlensing planet searches: high-magnification 
events that are intensively monitored over the peak. The observations of these events are 
always frenetic, sometimes even comical, so it may seem surprising that they nevertheless 
constitute a "controlled experiment", or very nearly so, and hence are subject to rigorous 
analysis. 



More than a decade ago lGriest fc Safizadehl ( 1l998l ) pointed out that high-mag events are 



much more sensitive to planets than their far more common cousins, the low-magnification 
events. The reason is fairly simple. For point-lens/point-source microlensing, the magnifi- 
cation is very nearly A — )■ for u < 1/3, where u is the source-lens separation normalized 
to the Einstein radius. "High-magnification" (Amax ^ 1) therefore implies impact param- 
eter uq = Aj^ax ^ 1' th^t source is projected very close to the lensing star. If the 
lensing star has a planet, the planet's gravitational field will perturb the otherwise circularly 
symmetric gravitational field of the star, and so induce a tiny "central caustic" (contour of 
infinite magnification) around it. If the source passes over or near this caustic, the lightcurve 
is perturbed, reveahng the planet's presence. The key point is that the planet can be at a 
broad range of separations and at virtually any angle relative to the source trajectory and 
still create a caustic that will perturb the light curve. 

However, following this seminal paper, microlensing follow-up monitoring groups con- 
tinued to focus primarily on garden variety microlensing events, which are only sensitive to 
so-called "planetary caustics". These are formed by the action of the primary-star grav- 
itational field on the planet gravitational field, and as such are bigger and so have larger 
cross sections than central caustics. At first sight, this makes them more favorable targets, 
but this conclusion only holds if one has unlimited telescope time for monitoring. Then 
one would monitor all available events (and hence primarily low-mag events) and would find 
most planets in these events (just because the caustics have a larger cross section). But 
if observing resources are limited, then one should focus these on high-mag events because 
these can be predicted (at least in principle) from the pre-peak part of the light curve and 
have individually higher sensitivity to planets. By contrast, there is no way to predict that 
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a source is approaching a planetary caustic. 



Nevertheless, within the context of continued focus on planetary caustics in normal 
events, there were significant efforts to take advantage of high-mag events as well. In 1998 
MPS issued the first high-mag alert for MACHO-98-BLG-35, which receive d an enthusiastic 
response from the MOA group, leading to the possible detection of a planet (IRhie et al.lbOOd : 
Bond et al.ll2002bl ). but with too-low significance to be confident. From the inception of its 
survey, MOA made real-time alerts of high magnification events a priority an d attempted 



to org anize follow-up from other continents, with 10 such alerts the first year (iBond et al. 



2002aD. The most spectacular success of this program was MOA-2003-BLG-32/OGLE-2003- 
BLG-219, which was densely sampled over its A^ax = 520 peak by the Wise observatory in 
Israel after such an alert, and which yielded the best upper limits on planetary companions to 
a lens to that date ( lAbe et al.ll2004j ). Theo retical work was also do ne to optimize observations 
of high-mag events for planet sensitivity f Rattenburv et aPbood ). 



When the Microlensing Follow Up Network (//FUN) began operations in 2001, it followed 
the already-esta blished model of follow-up observations, which did include high-mag events 
( lYoo et al.ll2003l ). but did not emphasize them. However, three things happened to change its 
orientation toward concentrating on high-mag events. First, the OGLE-III survey came on 
line in 2002 with a discovery rate of 350 events per year, moving up to 600 events per year in 
2004. This compared with 40-80 events per year discovered by OGLE-II in 1998-2000. The 
number of events alerted per year has a direct impact on whether or not one is in the regime 
of "limited" or "unlimited" resources. Before 2002, if one restricted oneself to high-mag 
eve nts, one would most ly be sitting on one's hands. For example, when lAlbrow et al.l (120011 ) 
and lGaudi et al.l (120021 ) analyzed 5 years of PLANET collaboration data, they reported only 
two events with Amax > 100. OGLE-III dramatically changed that situation. More recently, 
the MOA collaboration inaugurated MOA-II (2007 first full season), which has had the net 
effect of increasing the total rate of reported events by about 50%. 

Second, /xFUN began attracting the intrinsically "limited" observing resources of ama- 
teur astronomers. These contrast with the larger dedicated professional observatories in two 
key ways. First, the observers generally cannot observe all night, every night, or they will 
be unable to keep their day jobs. Second, the smaller apertures of their telescopes restrict 
them to relatively brighter targets. Both "limitations" naturally drive amateurs to high-mag 
events, which have a bigger chance of science payoff and are brighter (because highly mag- 
nified). Moreover, there is one crucial dimension in which amateurs are not limited: they 
have completely free and almost instantaneous access to their telescopes at any time. Thus, 
while dedicated follow-up telescopes are typically operated only in the 3-4 month core of the 
season (when microlensing targets are observable for at least half the night), amateurs can 
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react to alerts deep into the wings of the season, close to doubling the number of high-mag 
events that can be monitored. In 2004, one amateur began observing /iFUN targets on her 
own initiative. She requested regular alerts and began organizing other amateurs to join in, 
who in turn self-organized a network. About half of the /iFUN authors of this paper are 
amateurs. 

Third, /xFUN had to become aware of this changed situation, i.e., that it had moved from 
the domain of unlimited to limited resources. This transition was partly aided by the fact that 
yuFUN access to two of its professional telescopes (Wise and SMARTS CTIO) was limited 
by their being shared resources. But throughout 2003-2004, /xFUN "straddled two horses", 
focusing on high-mag events when available, but trying to keep to the old planetary-caustic 
strategy most of the time. Pr eparations for the 2005 season were significantly influenced by 



preliminary work (ultimately, iDong et al.l 120061 ) showing that unless high- mag events were 



intensively monitored over their peak, much of their sensitivity to planets is compromised. 

Then in April 2005, /iFUN intensively followed the (by those days' standards) high- 
mag ev ent OGLE-2005-BL G-071, which resulted in the detection of the second microlensing 
planet ( lUdalski et al ] l2005h . This detection led to /iFUN consciously changing its orientation. 



procedures, recruitment, etc., with the aim of focusing primarily on high-mag events. 

From 2005 onwards, considerable effort has gone into identifying potential high-mag 
events, and in some cases obtaining additional data to improve the prediction of Amax- If an 
event is deemed a plausible high-mag candidate, then observers are notified by email, without 
necessarily being urged to observe, just to put them on alert. Once high magnification 
seems probable, observations are requested with various degrees of urgency. The urgency is 
conditioned by the fact that peak sensitivity is usually less than a day (the normal human 
cycle time) but more than a few hours (so requiring observations from multiple continents). 
Many factors enter into the quality of the final light curve, including weather conditions on 6 
continents plus Pacific Islands, observer availability, communication problems, etc. Indeed, 
it is difficult to convey the level of chaos during one of these events. 

Despite (and also because of) this chaos, the resulting data stream generally retains the 
character of a "controlled experiment". In the "ideal incarnation" of this search mode, the 
event is recognized (with greater or lesser certainty) to be approaching a high-magnification 
peak, and an alert is issued to interested observers urging them to observe it intensively over 
the predicted peak. The observations take place regardless of whether the planet is present or 
not. The very chaos, remoteness of observing locations and communication problems make 
it difficult to gain knowledge of planetary perturbations until after the key observations are 
over. In Section [21 we discuss how well the real events conform to this ideal, but assuming 
for the moment that they do, the ensemble of high-magnification events with and without 
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detections can be analyzed in exactly the same way as can be done for the "pure survey" 
mode. In Section [31 we derive an analytic expression for the sensitivity of an ensemble of 
high magnification events monitored densely over their peak, and in Section HJ we calculate 
numerical sensitivities and use these to measure the frequency of ice giants and gas giants 
beyond the snow line. Finally, in Section [5l we discuss some implications of our results. 

2. Selection Criteria and Data 

We begin by designing criteria for selecting events and planets to be included in the 
analysis that enable the detections and non-detections to be analyzed on the same footing. 
This is essential to the goal of defining a data sample that can be treated as a "controlled 
experiment" . 

2.1. Selection Criteria for Events 

El) High-cadence (< lOmin) /iFUN data over some portion of the peak 

E2) High-cadence data covering at least one wing of the peak, |t — to| < ^eff, with no major 

gaps 

E3) High signal-to-noise ratio (S/N) data (Xli > 50000) covering the other wing, where 
the ai are the data-point errors in magnitudes 

E4) Light-curve is not dominated by binary (non-planetary) features 

We define tes = mo^e for point-lens/point-source events, where is the Einstein crossing 
time. For point-lens events that suffer finite-source effects, we generalize this definition to 
the time interval during which the magnification is within a/2 of the peak. And for planetary 
events, we further generalize it to the time interval that the magnification would have been 
within a/2 of the peak if the lensing star had lacked planets. 

We now justify these criteria. 

In addition to reflecting the fact that we are summarizing /xFUN work, criterion (El) 
ensures that we can rigorously review the available data on ~ 3000 events discovered during 
2005-2008, and reduce them to a manageable subset of "only" 315 that can be investigated 
using /iFUN files. These 315 were then quickly pared down to a few dozen events that are 
consistent with high-magnification and actually meet criterion (El). 

Criteria (E2) and (E3) should be considered together. Our underlying requirement 
is to have enough coverage of the event so that planetary deviations that give rise to a 
Ax^ = 500 deviation (for the rereduced data set), have a high probability of yielding a 



- 14 - 



unique scientific interpretation. See criterion (P2), below. These two c riteria are basicall y 



derived from the experience analyzing the event OGLE-2005-BLG-169 (IGould et al.ll2006l ). 
which does somewhat better than barely satisfy both. The analysis of this event was already 
fairly difficult because of the multiple minima and would have been quite degenerate if, 
for example, it lacked high S/N data on the rising side. Hence it would not have led to a 
publishable planet with a reasonably well-defined mass ratio. 

Criterion (E4) is adopted for two reasons. First, in contrast to point-lens events, most 
binary events are not modeled with sufficient precision to measure the uq parameter well 
enough to construct a well-defined sample of events with maximum magnification greater 
than some threshold Amax- Second, the problem of detecting planets in the presence of light- 
curve features dominated by a binary is not well-understood. Hence, the results we derive 
here really apply to stars not giving rise to strorig bin ary features, which excludes roughly 



3-6% of all stars flAlcock et all 120001 : 1.Taroszvhskill200ffl . 



Table [T] lists all events from 2005-2008 that satisfy these four criteria and that had 
peak magnifications Amax > 100. Column 1 is the event name, column 2 is the maximum 
magnification, column 3 is the time of closest approach between the source and lens to, 
column 4 is the Einstein timescale ^e, column 5 gives the mass of the lens star for cases that 
it is known, and column 6 gives the method by which it is derived. These lens masses and 
methods are discussed in Section 14. 2[ The events are listed in inverse order of Amax- 



Figure [H shows the cumulative distribution A^^^ for the 15 events in Table [H It displays 
a clear break at Amax = 200. Below this value, the distribution is uniform in ^^ix) which 
is the expected behavior for a complete sample (ignoring finite source effects), i.e., uniform 
in uq, which is the impact parameter in units of the Einstein radius. T he dashed lin e , with 



slope of (iiVev/(iA~ax — 2600 is a good match to the A^ax > 200 data. I Cohen et al.l (120 lOl ) 
found such uniformity for the underlying sample of OGLE events in 2008, with a slope of 
dNev/duo = 1080 for mq < 0.05. From this comparison, we learn two things. First, /iFUN was 
aggressive enough to achieve a uniform subsample only for events with Amax > 200. Second 
^FUN was able to intensively monitor half of all events in the A^ax > 200 subsample. 
That is, assuming that OGLE found similar numbers of high-mag events in 2005-2007, and 
accounting for the fact that MOA found 50% more events (not found by OGLE) in 2007- 
2008, the full sample of high-mag events was about dN^v/duo ~ 5 x 1080 = 5400, of which 
yuFUN effectively monitored about 48%. 

Although it may not be immediately obvious. Figure [1] implies that we must impose a 
fifth criterion. 
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E5) A^,, > 200 

Figure [T] demonstrates that /iFUN was substantially less enthusiastic about events 
^max < 200 than A^^x > 200, whether because it simply did not act on events known 
in advance to be in the former category or just became less enthusiastic about observations 
once these events were recognized near peak not to be extremely magnified. This bias is 
a natural consequence of /xFUN's limited observing resources (as discussed in Section [T]): 
there are 4 times as many events with Amax > 50 as A^^x > 200 and their duration of peak, 
2 teff ~ 2 tE^mL lasts 4 times as long, so 16 times more observing resources would be required 
to follow them all. Hence, if an event proved midway to be one that the objective evidence 
demonstrates /iFUN cared less about, there would be a tendency on the part of observers 
to slacken efforts (whether or not the internal alert was officially called off). Then the event 
would have less chance of meeting the selection criteria. But if a planet were detected during 
the peak observations (and there is a greater chance of recognizing a planet in real time for 
lower Amax because the peak lasts longer) then observations would not slacken, but rather 
intensify. Since this bias cannot be rigorously quantified, planets and non-detections from 
^max < 200 events must both be excluded from the sample. 

2.2. Selection Criteria for Planets 



PI) Planet must be discovered in an event that satisfies (E1)-(E5) 
P2) Planetary fit yields improvement Ax^ > 500 
P3) Planet-star mass ratio q must lie in the range 



10 



-4.5 



10" 



Criterion (PI) is self-evident but is stated explicitly for completeness and emphasis. 
Criterion (P2) may appear at first sight somewhat draconian, but it is realistic. To explain 
this, we first note that among the six planets liste d in in Table El the "weakest" detection 
is Ax^ = 880, which is for MOA-2008-BLG-310Lb JPong et al.lboOQbl l Now, it is certainly 
possible to recognize systematic residuals frora a poi nt-lens fits "by eye" at a much lower 
level, even Ax^ = 100. Indeed, iBatista et al.l (120091 ) argued that no systematic residuals 
were present in the fit to OGLE=2007-BLG-050 at a much lower level, Ax^ = 60. But 
if such deviations had been observed in an event, this would not have necessarily enabled 
discovery of a planet, where "discovery" here means "publication". First, Ax^ = 100 in the 
final, rereduced and carefully cleaned data implies something like Ax^ = 50 in the standard 
pipeline data, and systematic deviations due to a planet at this level would probably not be 
recognized as significant, i.e., clearly distinguishable from systematics that appear in dozens 
of other events, and which just refiect instrumental, weather, or data-reduction problems. 
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But more to the point, even if the unrereduced-data were Ax^ = 100, triggering strenuous 
efforts to clean and rereduce the dataset, resulting in, say, a Ax^ = 200 improvement, it 
is far from clear that this deviation (even if strongly believed to be real) would lead to a 
publishable planet detection. This is because, in addition to obtaining an acceptable fit to 
the data, such a paper would have to demonstrate that there could be no acceptable fits 
to the data to non-planetary solutions. We have already designed criteria (E2) and (E3) 
to eliminate those events for which very high Ax^ is possible without leading to a unique 
interpretation (due to incomplete coverage of the deviation). But we still must set the 
threshold high enough so that if an anomalous event survives criteria (E2), (E3), and (P2), 
it has a small chance of being ambiguous in its interpretation. Our best estimate of this, 
from experience fitting events, is Ax^ = 500. However, we regard other values in the range 
350-700 as also being plausible candidates for this threshold. We will show in Section H] that 
our basic conclusions are robust to changes within this range. 

The upper boundary in Equation ([1]), Criterion (P3), is necessary because at high mass 
ratios q, one cannot be confident that the event will not be rejected (consciously or un- 
consciously) as a "brown dwarf" or "low-mass-star", and therefore not be monitored as 
intensively as it might be [and so not pass criteria (E2) and (E3).] To illustrate this, we 
review how OGLE-2008-BLG-513Lb (Yee et al., in preparation) was "almost rejected" as a 
binary, even though it is probably a planet. This event has a large, strong, resonant caustic, 
that was initially mistaken for a binary. During the long intra-caustic period, it was realized 
that the companion might be a planet, and that intensive observations of the caustic exit 
would be necessary to resolve this question. Such observations were obtained, and from these 
we know that the impact parameter was uq ~ 0.027, so Amax ~ 37, which means that the 
event fails criterion (E5) and so is excluded from our sample. But if these data had not been 
obtained, then uq for this event would not be known, and it would not be known whether 
the event was in the sample or not, and if it were, whether the companion was a planet or 
not. 

Why does this example then not just prove that the whole concept of "controlled exper- 
iment" is unviable? The answer is given by Figures [2] and [31 below. One sees from these that 
at high mass ratios, q ~ 10~^, there is an extremely wide range of s for which the planet is 
"detectable". Only a small fraction of these "detectable" events have s ~ 1, which produce 
large resonant caustics that might be mistaken for binaries. Hence, while in principle some 
of these "detections" might be lost to this confusion, the great majority would not cause 
any confusion. The reason that planets like OB08513Lb make their way into the detections 
at all, despite their relative rarity, is that the caustics are so large that they are detectable 
over a wide range of uq, only a small fraction of which would pass the "high-mag" criteria 
of Section 12. 1[ Nevertheless, as q grows, this potential problem grows with it. We adopt 
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logg+ = —2, but recognize that values ranging from 
plausible choices. 



-2.3 to —1.8 might also have been 



The lower boundary g_ = 10~^'^ is established because of concerns of the real "de- 
tectability" of lo w mass planets in the presence of higher mass planets in the same system. 
In the method of iRhie et al.l ( 120001 ). which we employ in Section W7L\ planet sensitivity is de- 
termined by fitting simulated star-planet light curves that are constructed to have the same 
error properties as the actual data, to point-lens models. If the best such model increases 
by more than a given threshold (say Ax^ > 500), then the planet is said to be detectable. 
Since this method directly mimics the process of planet detection for single-planet systems, 
it is a good way to characterize the detectability o f such systenis. B ut high-magnification 
events are particularly sensitive to multiple planets (jCaudi et al. 1998 ) and why should this 
approach tell us anything about the det ectability of a seco nd planet in a system already 
containing one planet? As first shown by iBozza et al.l (119991 ). the net perturbations of such 
two-planet high-magnification lightcurves usually "factor" into the sum of perturbations in- 
duced by each planet separately. For example, th e only published two-planet system has this 
property ( iGaudi et al.ll2008l : iBennett et al.ll2010l ). Of course, the factoring is not perfect, but 
in this real case (and in many simulated cases), once the dominant-planet perturbation is 
removed, the secondary perturbation is easily recognized, leading to an excellent starting 
point for a combined fit to both planets simultaneously. For reasonably comparable planet 
mass rat ios, the only exception to this is i f the planet-star axes are closely (within < 20°) 



aligned (IRattenbury et al 



2002; Hanlboosi ). In this case, the single-planet fit still fails, but 



the residuals to this fit are not easily recognizable. While such difficulties might impede 
recognition of the second planet, the required alignment is so close, that such cases would 
be a small minority of two planet systems. 

However, t his factorin g has only been studied in detail for planets with relatively com- 
parable masses I2OO5I ) . The situation may not be as simple when the mass ratio of the 
two planets is extreme. Based on analysis of the events listed in Table O below, we cannot 
be confident of excluding all "second planets" with q < 10~^. To be "conservative", we have 
moved the boundary to g_ = 10~^'^. 

Because both g_ and g+ have some uncertainty, we must ask how robust our conclusions 
are to changes in these parameters within a reasonable range. We show in Section H] that 
our basic conclusions about planet frequency are not seriously affected by uncertainty in q_ 
and g+. However, these uncertainties will prevent us from deriving a slope of the mass-ratio 
function. 

Six planets satisfy criteria (P1)-(P3). Table |2] displays their characteristics. Columns 
2-5 are the parameters {uq, p,q,\\ogs\) measured from the event, where p is the source 



- 18 - 



radius in units of the angular Einstein radius 6^. In several cases, there is an unresolved 
s ^ s^^ degeneracy, which is irrelevant to the current study, so we just display the absolute 
value of the log. The final column is the maximum detectable value of | logs| according to 
calculations reported in Section 14. 1[ 



2.3. Were Discovery Observation Cadences Really Independent of the Planet? 

Were the observations that led to the discovery of the six planets in fact carried out 
i ndependent of th e presence of the planet? For three of the se, OGLE-2005-BLG-169Lb 
(IGouM et allbcjOeh. M OA-2007-BLG-400Lb flOong et allboOObh . and MOA-2008-BLG-310Lb 



( iJanczak et al.ll2010l ). the planet was not recognized until after the event had returned to 
baseline. OGLE-2007-BLG-349 (Dong et al., in preparation) was recognized to have a sig- 
nificant deviation possibly due to a planet based on observations in Chile, 36 hours after the 
call for intensive observations based on its high-mag trajectory, and roughly 7 hours after 
observations had begun in South Africa. While it is true that reports of this potential planet 
heightened excitement, and may perhaps have increased the commitment of observers to get 
observations, the density of observations (from 4 continents plus Oceania) did not qualita- 
tively change from before till after the potentially-planetary anomaly was recognized. This is 
the only one of the six planets in our sample that is not yet published. The reason is that the 
system contains a third body, which has proven difficult to fully characterize. However, the 
characteristics of OGLE-2007-BLG-349Lb are very well established, and the third body is 
certainly not in the mass range being probed in the current analysis. Hence we feel confident 
including this planet in the sample. 



The two planets OGLE-2006-BLG-109Lb,c f lGaudi et al.ll2008l : iBennett et al.ll201(l ) re- 
quire closer examination. OGLE-2006-BLG-109 was recognized to be an interesting event 
almost 10 days prior to peak due to detection of what turned out to be the resonant caustic 
of OGLE-2006-BLG-109Lc, the Saturn mass-ratio planet in this system. This anomaly did 
indeed trigger some additional observations, which did help characterize this planet. But 
a review of email communications that initiated follow-up observations during the event 
shows that far more intensive observations were triggered several days later, after the event 
had appeared to return to "normal" (point-lens-like) microlensing, exactly by its high-mag 
trajectory. Although these emails remark on the possible presence of a planet, they place 
primary emphasis on this being an otherwise "normal" microlensing event that was reaching 
extreme magnification. (Note that the appearance of extreme magnification was not itself 
an artifact of the presence of planets(s), but was simply due to low source-lens impact pa- 
rameter.) It was the intensive observations from New Zealand, triggered by these emails, 
that captured the "central structure" of the caustic due to OGLE-2006-BLG-109c. These 
would have enabled basic characterization of this planet even without the fiurry of followup 
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observations 10 days earlier. Moreover, it was the same email that triggered intensive ob- 
servations from two widely separated locations (Israel and Chile), that enabled detection of 
the "central caustic" due to OGLE-2006-BLG-109Lb, the Jupiter mass-ratio planet. Thus, 
all detections were in reasonable accord with the "controlled experiment" ideal. 

3. Analytic Treatment 
3.1. Triangle Diagrams 



Batista et al.l (120091 ) and lYee et al.l (120091 ) recently analyzed the sensitivity to planets 
of two events in our sample from Table [H OGLE-2007-BLG-050 and OGLE-2008-BLG-279 
respectively. Figures [2] and [H] are versions of their results, but with the detection thresh- 
old used in this paper, Ax^ = 500. In contrast to the wide range of detection-sensitivity 
morphologies shown in Figure 8 iGaudi et al.l (l2002l ) . both of these diagrams have a simple 
triangular appearance, which is basically described by a two-parameter equation. 



log S I max (g) = r/log 



(2) 



where q is the planet/star mass ratio, s is the planet-star projected separation in units of the 
Einstein radius, r] is the slope of the triangle, and q^i^ defines the "bottom" of the triangle. 
Planets lying inside the triangle are detected with 100% efficiency and those lying outside 
are undetectable. The boundary region is quite narrow. In principle, planet detection is 
a function not only of (s,g), but also a, the angle of the source-lens trajectory relative 
to the planet-star axi s. The narrowness o f the boundary reflects that detection is almost 
independent of a. See iBatista et al.l (120091 ) and lYee et al.l (120091 ) for concrete illustrations. 



It is also striking that the slopes t] ~ 0.32 and rj ~ 0.35 are nearly the same for the 
two diagrams, leading to the conjecture that rj is very nearly co nstant for wel l -moni tored 
high-magnifi c ation events. Indeed, of the 43 events analyzed by lAlbrow et al.l ( 120011 ) and 
Gaudi et al.l (120021 ) . two are relatively high-mag (Amax > 100) and both have the same tri- 
angular appearance and very similar slope rj, as does t he extreme -A max = 3000 event OGLE- 



BLG-2004-343 with simulated coverage analyzed by iDong et al.l (120061 ). If truly generic, 
this would mean that high-mag event sensitivities have an extremely simple triangular form 
characterized by a single parameter, gmin- 

Moreover, while gmin, i-e. the depth of the triangles in Figures |2] and [3l obviously depends 
on the intensity, quality, and uniformity of coverage, one expects the fundamental scaling to 
be 

Q'min = ^^max' (3) 

where ^ is a parameter that depends on the data quality, etc. The reason for this expected 
scahng is that the size of the central caustic is proportional to q, and v4~^ measures how 
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closely the source probes the center, which i s roughly the maxim um of the impact parameter 
Uq and the source size p ( Han fc Kim 2009 : Batista et aLll2009l ) 



Hence, armed with an empirical estimate of ^, one can quickly gage the sensitivity of one 
event or an ensemble of events to planets, which is quite useful both to guide observations and 
as a check on "black-box" simulations of event sensitivity. Indeed as we will show below, one 
can approximately "read off" the frequency of planets by just counting the number detected 
and the number of high-mag events surveyed. 

Based on this handful of published analyses, we estimated C, ~ 1/70 for a Ax^ = 
500 threshold. Since these analyses (naturally) focused on events with better-than-average 
coverage, we estimate that C, ~ 1/50 is more appropriate for a sample such as ours. 



3.2. Analytic Estimate From Triangles 



We now assume that planets are distributed uniformly in logs [Opik's law f lOpiklll924l )] 
in the neighborhood of the Einstein rin^. We will show further below that this assumption 
is consistent with current microlensing data. Then, assuming rj = 0.32 to be universal, we 
have 

q 

Pi{q)d\ogq = 2r]fi{q) log e(g - qmin,i)d\ogq (4) 

Q'min,i 

where /(g) is the number of planets per dex of projected separation per dex of mass ratio, 
and 6 is a step function. The expected number of planets detected in high-mag events is 
then just the sum of Equation (jl]) over all high-mag events with good coverage 



dK ^' 



'pi 



dlogq 



Y.P^{q)- (5) 



i=l 

Figure [1] demonstrates that the /iFUN sample is uniform for < A~^^-^ < e, where e = 0.005. 
Hence, we can turn this sum into an integral. 



dlogq 



Y,P^{q)^— / ^A_(g)rfA^L = 2r//(g)— / log[g/(e/A^a.)]e(g-e/A^ax)])rfA-L 



f{(l)— log(g/gmm)0('? - gmin)c?gmm = 77T/(g)— / \n{q / q^i^)dq^i^, 

e Jo 4 In 10 e Jq 

(6) 



^ If in fact planets are distributed dN/d\ogs oc s^, then Eq. ^ is in error by sinhx/x with x — 
?7plng/(7niin. For ry — 0.32, p — 0.4 and q/qmin — 100, this is still only a factor 1.06. 
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which may be evaluated, 

ctlogg mlO 

9{q) = — (g<ee); (?(g) = l + ln^ (g>ee) (8) 

qthv Qthr 

where gthr = ^e- Below this threshold, detection efficiency falls linearly with q. Above 
the threshold, it rises logarithmically with q. Note that the appearance of "In 10" in these 
formulae is an artifact of our having chosen to express the density of planets in units of dex 
of separation, rather than the "natural unit" of an e-folding. 

Both the break at gthr in the normalized survey sensitivity g{q) and the functional forms 
of g{q) on either side of this break are easily understood from the triangular form of the 
individual-event sensitivity diagrams. For q > q^hr, all N^v events contribute sensitivity. If 
we compare two mass ratios, logg and logg + rflogg, the latter is sensitive to a log-separation 
interval on the triangle that is larger by exactly 2rid\ogq for each individual event, so the 
sensitivity of the ensemble of events is simply 2?7 log g+const. On the other hand, for q < gthr, 
only a fraction g/gthr of the events contribute, which breaks the logarithmic form of g{q)- 

Note that our estimate ^ = 1/50 implies that our Amax < = 200 survey has 

gthr = = 10-', (9) 

i.e., twice the mass ratio of Neptune. We discuss the implications of this threshold in 
Section |5l 



4. Numerical Evaluation 



4.1. Sensitivities for the Full Sample 



To more accurately determine the sensitivity of our survey and to infer the frequency 
of planets, we carry out a detailed sensitivity analysis of all 13 of the A^ ^^ > 200 ev e nts in 



Table [H (except the two that were already done). We use the method of iRhie et al.l f|2000l') 



(outl ined in Section 12. 2j) except that we take full account of finite source effects (IDong et al. 
20061 ). which are much more important for the present sample of events because of their 



higher magnification. For the last three events above the cut in Tabl e [H 6p, (and hence 
p = 9^/9e) is not well measured. For these, we follow the procedure of iGaudi et al.l (120021 ) 
and adopt p = 6^./{ptypt'E), w here 6^, is deterrn ined from the instrumental color- magnitude 
diagram in the standard way f lYoo et al.ll2004l ). is the measured Einstein timescale, and 
/ityp = 4masyr~^ is the typical source-lens proper motion toward the Galactic bulge. Figure 
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m is a "portrait album" of the resulting triangle diagrams (only one side shown to conserve 
space) and Figure E] shows the integrated sensitivity of each event as a function of mass ratio 
q. That is, it is the integral of the sensitivity (in Figs. E], |3], and H]) over horizontal slices. 
Hence, if the sensitivity were truly a triangle, the curves in Figure [5] would be perfectly 
straight lines, with slope 2ri (illustrated by the bold black hne segment) and x-intercepts at 
gmin- Most of the curves do have this behavior over the range —4 < logg < —2, which is the 
main range of sensitivity of this technique and also where the planets in Tableware located. 
Moreover, the inferred intercepts of the straight-line portion of these curves do generally 
reach to lower mass ratio gmin for higher Amax events, although with considerable scatter. 
However, while some events (like OGLE-2008-BLG-279) are almost perfectly straight down 
to zero, others (like OGLE-2005-BLG-169 and OGLE-2006-BLG-109) show a sharp flattening 
toward lower mass ratios. There are two reasons for this. Events like OGLE-2005-BLG-169 
have non-uniform coverage over peak, which makes detectability a strong function of angle. 
The contours in the triangle diagram separate, so that while the 50%-sensitivity contour is 
fairly straight, there is still substantial sensitivity below qmm, which is defined by where the 
two 50% contours meet, creating a long tail of sensitivity below this threshold. Events like 
OGLE-2006-BLG-109 have very small source size relative to impact parameter, p/uq <^ 1, 
which enables detection of small mass-ratio planets that are very close to Einstein ring 
because the relatively large, but ve ry weak, caustics of th ese planets are then not "washed 
out" as they would be for larger p teennett fc Rhiel l Il996l ). Hence the entire "triangle" has 



a curved appearance, although the contours are tightly packed together. 

The black bold dashed curve in Figure [5] is the combined sensitivity, i.e., the sum of the 
sensitivities for all 13 events (divided by 10, so it fits on the same plot), which we call G{q). 

The curves in Figure allow us to compare the observed log projected separation s, 
with the maximum detectable separation. See Table IH Because some planets suffer from the 
s ^ s^^ degeneracy, we only show the absolute value of log s. The cumulative distribution 
of the ratios of these quantities is shown in Figure [61 The separations are consistent with 
being uniform in log s, with Kolmogorov-Smirnov probability of 20%. (Moreover, in general, 
the high-magnification events with the greatest values of | logs| - and so the greatest po- 
tential leverage for probing the distribution as a function of s - are also the most severely 
affected by the s degeneracy. This applies, for example, to MOA-2007-BLG-400Lb. 

Thus the dependence on s will be much better explored using "planetary c austics" in low 



magnification events, for w hich the s ^ s ^ degeneracy is easily resolved. iGould fc Loeb 
I992I : iGaudi fc GouldlEgg?! ) 



4.2. Masses of Host Stars 



Now, f{q) may in principle be a function of the host mass M (and perhaps other 
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variables as well). With only six detections, we are obviously in no position to subdivide 
our sample. Nevertheless, it is important to assess what host mass range we are actually 
probing. There do exist mass estimates or limits for all 5 hosts of the planets that have 
been detected, and there are also mass estimates for 5 of the 8 lensing stars in the sample 
for which no planet was detected, which are given in Table [1] together with the method 
of estimation. For 5 of the events, both the angular Einstein radius 6e and the "microlens 
parallax" tte are measured, whic h together per mits a mass measurement M = 6-^/k,7t-e, where 



K = 4G/(c2AU) ~ 8.1 mas Mq^ flGouldlbonnbl ). The measurement for OGLE-2007-BLG-349 



is preliminary (Dong et al., in preparation) but the others are secure. For four of the events, 
there is a measurement of 6e = ^y nMTTrei, where vTrei is the source-lens relative parallax, 
but not tte. This measurement of the product of M and TTi-d combined with the lens-source 
relative proper motion fi = 6-E/t-E and a Galactic model (GM) permit a Bayesian estimate 
of M. Finally there are two events for which adaptive optics (AO) observations provide 
information on the host. For OGLE-2006-BLG-109, AO resolution of the host confirms the 
microlensing mass determination from 6e and tte. For MOA-2008-BLG-310, AO observations 
detect excess light (not due to the source) but it is not known whether this excess is due to 
the lens or another star. So only an upper limit on the mass is obtained. 

Seven of these 10 measurements are in the range of middle M to middle K stars, while 
one is a brown dwarf and another is likely to be a late M dwarf. They cover a fairly broad 
range approximately centered on 0.5 Mq, with a tail toward lower mass. Of course, there are 
also 3 lenses in the A^ax > 200 sample for which there is no mass measurement or estimate. 
These have timescales of 10, 26, and 59 days, which is much more typical of microlensing 
events than the sample having mass measurements (with 3 events having tE > 100 days). If 
these are otherwise typical events, then the lenses probably lie mostly in the Galactic bulg e 



(IKiraga fc Paczyhskilll994j ). in which case their mean mass is roughly 0.4 Mq (lGouldll2000al ). 
Given that this is a minority of the sample, that the information about this minority is far 
less secure, and that the difference from the sample with harder information is not very large, 
we adopt 

M~O.5M0 (10) 

for the typical mass of the sample. However, we note that the implications discussed in 
Section |5] would not be greatly affected if we had adopted M ~ 0.4 Mq. 

4.3. Likelihood Analysis 

To evaluate the mass-ratio distribution function /(g) = Ag", we maximize the likelihood: 

Iv = -A^exp + X!lnG(gi)/(gO; iVexp = H dqG{q)f{q) (11) 
i=i -^1- 
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and find 



ct log get logs \5 X 10 



However, as we now argue, while the normalization of Equation (IT2l) is robust, the slope is 
not. 

In Section 12.2^ we summarized why there must be some boundaries q± beyond which 
the experiment is seriously degraded, but argued that there is no "impartial algorithm" for 
deciding exactly where those boundaries should be. We find that if we vary g+ between 
—2.3 and —2., and we vary g_ between —5 and —4.5, that the normalization in Equation 
(!T2|) varies by only ±10%, which is much smaller than the statistical errors. However, the 
power-law index varies between —0.6 and —0.2. If we had a much larger sample of planets, 
then we could set the boundaries at various places within the range of our detections, thereby 
simultaneously reducing both the number of detections and A^^exp in Equation ( ITTl) . For an 
infinite sample, such a procedure should lead to no variation in either slope or normalization. 
For a finite sample, the variation would provide an estimate of the error in these quantities 
due to the uncertainty in knowledge of these boundaries. However, when we apply this 
procedure to our small sample, we find quite wild variations, implying that we cannot derive 
a reliable slope from these data. 

In Section 1212] we mentioned that the threshold value Ay^ > 500 also had some intrinsic 
uncertainty. However, in this case the effect is very small. For example, if we decrease the 
threshold to Ax^ > 350, then the normalization in Equation (fT2ll decreases by only 7%, 
much less than the Poisson error. Given that the normalization in Equation f|T2l) is robust 
but the slope is not, we give our final result as: 

= (0.36 ±0.15) dex"^ at g ~ 5 x 10"^ (13) 

a log q a log s 

Figure [7] summarizes the principal inputs to the modeling. The top panel shows the 
cumulative distribution of the detections. The bottom panel shows the sensitivity of the 
survey as function of planet mass, both the analytic approximation derived in Section 13.21 
and the numerical determination derived in Section 14.11 These hardly differ. 



The density given in Equation ( 1131) . can be obtained by a very simple argument. The 
"triangle" for each event has sensitivity to (1/2 base x height) = ?7[log(0.01 Amax/O]^ ~ 
1.7[1 ±0.42 log(74max/400)]^ dex^ of planet parameter space. We observed = 13 events and 
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found six planets, so 6/(13 x 1.7) ~ 0.27 dex ^, i.e. correct within the statistical error. 

5. Discussion 



We have presented the first measurement of the absolute frequency of planets beyond 
the snow line over the mass-ratio range —4.5 < logg < —2. The resulting planet frequency, 
Equation ( lT3l) . can be understood directly from the data and the "triangle" sensitivity dia- 
grams. The result is applicable to a range of host masses centered near M ~ 0.5 Mq. The 
distribution is consistent with being flat in log projected separation s. 

5.1. Comparison with Previous Microlensing Results 



Gould et al.l ( l2006l ) had earlier concluded that "cool Neptunes are common" based on 



one of the planets analyzed here (O GLE-2005-BLG-169L b) and another planet with similar 
mass ratio, OGLE-2005-BLG-390Lb jBeaulieu et al.lbood ). which had been detected through 
another channel: foUowup observations of low magnification events. OGLE-2005-BLG-390Lb 
was actually recognized as a possible planetary event during the planetary deviation, but 
detailed review of these communications and their impact on the observing schedule shows 
that this "feedback" was not critical to robust detection. Moreover, at that time there 
had been no other detections through this channel, so the iBeaulieu et al.l ( l2006l ) detection 
could also be treated as a "controlled experiment" . These could then be combined to obtain 
an absolute rate for "c o ol Ne ptunes", albeit with large errors. However, with only two 



detections, iGould et al.l (120061 ) were not able to specify the mass range of "cool Neptunes" 
and hence were not able to express their results in units of dex~^ as we h ave done here. If 



we ne vertheless, somewhat arbitrarily, say that the iBeaulieu et al.l (120061 ) and iGould et al. 
(120061 ) result applies to 1 dex in logg, centered at the g = 8 x 10~^ of their two detections 



and adopt their 0.4 dex interval in log s, their estimated density of cool Neptunes can be 
translated to 0.95^0:55 dex ^ (lo")- To make a fair comparison with Equation ( |T3|) . it is 
necessary to adopt some slope for the mass function in order to account for the factor ~ 6 
difference in mass. If we adopt the lSumi et al.l (120101 ) slope of n = —0.68 fr om microlensing. 



then o ur prediction for this mass range would be 1.25 dex ^. If we adopt the lCumming et al. 



mm slope of n = -0.31 from RV, it would be 0.64dex-l Either way, these are consistent. 



Sumi et al.l (I2OIOI ) analyzed all 10 published microlensing planets, including the 5 that 



we analyze here. They approximated the sensitivity functions of this heterogeneous sample 
by a single power law (oc g*", m = 0.6 ±0.1) and derived a power-law mass-ratio distribution 
dNpi/dlogq q", n = —0.68 ± 0.2. Since we are unable to derive a slope from our analysis. 
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and they do not derive a normalization, there can be no direct comparison of results. 
5.2. Comparison with Radial Velocity Results 



Based on an analysis of RV planets, ICumming et al.l (120081 ) derive a normalization of 
0.035 dex~^, a factor 10 ± 4 smaller than the one found here. A factor (5 x 10~^/1.66 x 
10~^)~° '^^ = 1.5 of this difference is due to the fact that they normalize at higher mass 
ratio. The remaining factor 7 ± 3 difference is most likely due to t he different star-planet 



separations probed by current microlensing and RV experiments. The lCumming et al.l (12008 



study targets stars with periods of 2-2000 days, corresponding to a mean semi- m ajor- axis of 
a = 0^31 AU. Microlensing probes a factor ~ 3 beyond the snow line (Fig. 8 from lSumi et al. 



2010 



To make contact between microlensing observations of primarily lower mass stars 
with RV observations of typically solar type stars, we should consider planets in similar 
physical conditions, which we choose to normalize by the snow line. That is, we should 
compare to G-star planets at 3 "snow-line radii", i.e., a ~ 8 AU. Hence, the inferred slope 
between the RV and microlensing measurements is din N/d In a = log(7 ± 3)/log(8/0.31) = 
0.56 =b 0.16, wh i ch is consistent with the slope of din N/ din a = 0.39 ± 0.15 derived by 
Gumming et al.l (120081 ) for RV stars within their period range. Thus, simple extrapolation 
of the RV density profile derived from planets thought to have migrated large distances, 
adequately predicts the microlensing results based on planets beyond the snow line that are 
believed to have migrated much less. See Figure [HI Figure |9] compares microlensing and RV 
detections ClS db function of mass ratio q. 

5.3. Prospects for Sensitivity to Very Low Mass Planets 

Equation ([8]) and Figure [5] show a break in sensitivity at gthr — 10~^. For a power-law 
mass-ratio distribution, the ratio of planets expected above and below this threshold is 



Npijq > Qthr) 
Npiiq < Qthr) 



n + 1 



n 



2"ln(2) 



n 



n 



1) 



(14) 



where z = 0.01/gthr- This ratio is an extremely strong function of the adopted slope of 
the mass function, n. For n = —0.68 (ISumi et al.ll2010l ). it is 1.0, whereas for n = —0.31 



^ Just as RV measurements respond to a projected stellar velocities, and so measure msin i of the planet 
which is always less than or equal to the planet mass m, so microlensing observations measure the projected 
separation s, which for circular orbits is related to the semi-major axis by R^s = a sin 7 where 7 is the angle 
between the star-planet axis and the line of sight. The statistical distribution sin 7 is exactly the same as for 



sini in RV. Hence, except for rare cases when the orbit is constrained by higher-order effects (|Dong et al 



2009a: 



8 of 



Bennett et al. 2010f ). a must be statistically estimated from s (and Re), which is what is done in Fig. 



Sumi et al 



(201 
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(ICumming et al.ll2008l ). it is 4.7. The fraction of planets within the lower domain that lies 



below q is simply g""*" . Hence, while several authors have shown that i ndividual planets 



at or near Earth mass ratio are detectable in high-m agnification events (lAbe et al.l l2004l : 



Dong et al.l l2006l : lYee et al.l l2009l : iBatista et al.ll2009l ). the actual rate of detection will be 



strongly influenced by the actual value of n. As discussed in Section 12. 2^ probing to lower 
masses will require technical advances to robustly identify low mass planets in the presence 
of higher mass planets. But it will also require increasing the number of events that are 
monitored. 

There is some potential to do this. First, as shown in Section |U only about half of 
events that are announced by search teams are intensively monitored. Hence, there is room 
to double the rate by more aggressive monitoring. This will be aided by inauguration of 
OGLE-IV, which will have much higher time sampling and so will permit more accurate 
prediction of high-mag events. Second, it is possible that the more intensive OGLE-IV 
survey will increase the underlying sample of high-mag events. Finally, systematic analysis 
of high-mag events could bring down the effective Ax^ threshold from 500 to, say, 200, 
which would decrease ^ (and so gthr) by a factor (500/200)^/^ ~ 1.8. This would bring only a 
modest (logarithmic) increase in the sensitivity in the range q > qthr but would aid linearly 
for q < gthr- 

5.4. Constraints on Migration Scenarios 

We showed in Section |5^ that the planet density derived here, dNpi/d\ogqd\ogs = 
0.36 dex~^ is consistent with the density derived from RV studies, if the latter are extrap- 
olated to ~ 25 times the semi-major axis where the measurement is made. Regardless of 
the details of this comparison, the fact that the density of planets beyond the snow line is 
7 times higher than that at 0.3 AU, indicates that most giant planets do not migrate very 
far. Moreover, the fact that the slope found in RV studies at small s adequately predicts the 
density at large s, would seem to imply that whatever is governing the amount of migration 
is a continuous parameter. That is, it is not the case that there are two classes of planetary 
systems: those with migration and those without. Rather, all systems have migration, but 
by a continuously varying amount. This picture would be in accord with the evolving view 
of the Solar System, that even though the giant planets are in the general area of their birth 
"beyond the snow line", they have migrated to a modest degree. 

5.5. Comparison to Solar System 

Another interesting point of comparison is to the planet density in the Solar System, 
where there are 4 planets in the mass-separation regime that microlensing currently ex- 
plores, Jupiter, Saturn, Uranus, and Neptune. How common are "Solar-System analogs". 
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.e., systems with several giant planets out beyond the snow line? 



To address this question, we ask what the result of our study would have been if every 
microlensed star possessed a "scaled version" of our own Solar System in the following sense: 
4 planets with the same planet-to-star mass ratios and same ratios of semi-major axes as the 
outer Solar-System planets, but with the overall scale determined by the "snow line". While 
there is observational evidence from the ast eroid belt that the Solar-System snow line is 



near -Rsnow,© = 2.7 AU flMorbidelli et al.ll200Clf). there is considerable uncerta inty on how this 



scales with stellar mass (ISasselov Sz Lecarl 120001 : Kennedy &: Kenyoru l2008l ) . We therefore 
parametrize this relation by 

-^snow(-^^) -Rsnow,0 

hr^ h (15) 

and consider a range 0.5 < u < 2, adopting i/ = 1 for our fiducial value. We consider that 
the typical Einstein radius is Re = 3.5 AU (M/Mq)^/^ and the typical lensed star is M = 
0.5 Mq. Then scaling down the semi-major axis of a Jupiter analog so ajup-anaiog/-Rsnow = 
ajup/-Rsnow,0 implies 

'^Jup— analog '^Jup / \ 

" 3.5AUVMQy ' ^ ^ 

and similarly for the other 3 planets. We then imagine that these systems are viewed at 
random orientations, with the individual planets in random phases. From a Monte Carlo 
simulation, we find that we would then have expected (for our fiducial z/ = 1) to have detected 

18.2 planets (11.4 "Jupiters", 6.4 "Saturns", 0.3 "Uranuses", 0.2 "Neptunes") including 6.1 
systems with 2 or more planet detections. See Figure [101 For 1 < u < 2, the planet totals and 
multi-planet detections barely change. However they fall somewhat for smaller z/, reaching 

15.3 planets and 3.8 two-planet systems at u = 0.5. 

These results compare to 6 planets including 1 two-planet system in the actual sample. 

Hence, our Solar System appears to be 3 times richer in planets than other stars along 
the l ine of sight towar d the Galactic bulge. The single detection of a multi-planet sys- 
tem JCaudi et all l2008h allows the first estimate of the frequency of stars with "solar-like 



systems", defined as having multiple giants in the snow zone: 1/6, albeit with large errors. 
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Table 1: MONITORED Events with Magnification A > 100 



Name 




io(HJD) 


tE 


M/Mq 


Method 


OGLE-2007-BLG-224 


2424 


4233.7 


7 


0.056 ± 0.004 


M = 




OGLE-2008-BLG-279 


1600 


4617.3 


101 


0.64 ±0.10 


M = 


0e/ KTlE 


OGLE-2005-BLG-169 


800 


3491.9 


43 


0.49l°:i 


GMg 


© tE 


MOA-2007-BLG-400 


628 


4354.6 


14 


n on+0 19 


GMg 


© tE 


OGLE-2007-BLG-349 


525 


4348.6 


121 


~ 0.4 


M = 


QeI KTTe 


OGLE-2007-BLG-050 


432 


4222.0 


68 


0.50 ±0.14 


M = 


Oe/ ktte 


MOA-2008-BLG-310 


400 


4656.4 


11 


< 0.67 ±0.14 


AO 




OGLE-2006-BLG-109 


289 


3831.0 


127 


'-'•^-'--0.04 


M = 


^e/^tteAO 


OGLE-2005-BLG-188 


283 


3500.5 


14 


1 6+°-^^ 


GMg 


)^E ® ^E 


MOA-2008-BLG-311 


279 


4655.4 


18 


'-'■^'-'-0.09 


GMg 


)^E © ^E 


MOA-2008-BLG-105 


267 


4565.8 


10 








OGLE-2006-BLG-245 


217 


3885.1 


59 








OGLE-2006-BLG-265 


211 


3893.2 


26 








OGLE-2007-BLG-423 


157 


4320.3 


29 








OGLE-2005-BLG-417 


108 


3568.1 


23 









Tabic 2: Planets in Densely Monitored High-mag Events 



Name 


log (/() 


log 


log q 


1 log'H 


1 ( 1"' S 


OGLE-2005-BLG-169Lb 


-2.9 


-3.4 


-4.1 


0.009 


0.19 


OGLE-2006-BLG-109Lb 


-2.5 


-3.5 


-2.9 


0.20 


0.39 


OGLE-2006-BLG-109LC 


-2.5 


-3.5 


-3.3 


0.017 


0.27 


OGLE-2007-BLG-349Lb 


-2.7 


-3.3 


-3.5 


0.099 


0.48 


MOA-2007-BLG-400Lb 


-3.6 


-2.5 


-2.6 


0.47 


0.55 


MOA-2008-BLG-310Lb 


-2.5 


-2.3 


-3.5 


0.035 


0.14 
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Fig. 1. — Cumulative distribution of inverse maximum magnification for high-mag 

events observed by //FUN during the years 2005-2008. The distribution is uniform in A^^^ 
for high-mag events Amax > 200, with a slope of dNe^/dA^^ = 2600 {dashed line). A 
"controlled experiment" therefore requires a selection criterion Amax > 200. See Section [2?n 
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r — projected planet/star separation (AU) 
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Fig. 2.— Triangle diagram for OGLE-2007-BLG-050 (adapted from Fig. 9 of Batista et al. 



20091 ) showing sensitivity to planets of planet/star mass ratio q and planet-star projected 



separation d (alias s in the current paper). Planets within the contours would be detectable 
at Ax^ > 500 for 25%, 50%, 75% and 90% of source-lens trajectories. Since these trajectories 
are random, the contours reflect the probability of detecting the planet at a given s and q. 
The slope of the contours is t] = 0.32. 
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Fig. 3. — Triangle diagr am of planet sensitivities for OGLE-2008-BLG-279 (adapted from 
Fig. 7 of lYee et al.l[2009l l Similar to Fig. |3] except with contours of A^^ > 500 detectability 
with 10%, 25%, 50%, 75% 90% and 99% probability. Contour slope is = 0.35. 
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Fig. 4. — New "triangle diagram" sensitivities for all events in our sample except the two 
shown in Figures |2] and [3l Sensitivity is the fraction of all trajectory angles a that a planet 
would have been detected at Ax^ > 500, if a planet of mass ratio q had been present at 
projected separation s (in units of the Einstein radius). Only half of the diagram (which is 
almost perfectly symmetric) is shown here to conserve space. Positions (logg, | logs|) of all 
detected planets from Table |2] are shown as gold stars. 
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Fig. 5. — Calculated sensitives for each of the 13 microlensing events. These are integrals 
over horizontal cuts (at fixed q) in Figures |2], |3], and HI Events are shown in "rainbow order" 
according to magnification, with the first 7 in solid lines and the last 6 in dashed lines. The 
bold dashed black curve represents the total sensitivity of the sample (divided by 10). 
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Fig. 6. — Cumulative distribution of the (absolute value of the log of the) measured projected 
separation, | logs|, in units of the maximum value of this parameter, derived in Section |4?T] 
The observed distribution is consistent with planets being distributed uniformly in logs, 
with Kolmogorov-Smirnov probability of 20%. 
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Fig. 7. — Principal inputs to the modeling. Top panel: cumulative distribution of planet 
detections in 4 years of intensive monitoring of high-magnification events. Bottom panel: 
Sensitivity of the survey as a function of planet / star mass ratio q, both the analytic approx- 
imation derived in Section 13.21 and the numerical determination derived in Section 14.11 and 
shown in Figure |5l These hardly differ. 
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Fig. 



Planet f requencies determined from microlensing (this paper) and RV 



(jCumming et al.l |2008| ) at different semi-major axes. The RV result is scaled to mass ra- 
tio g = 5 x 10~^, using the RV-derived slope, n = —0.31. In order to take account of the 
different host star masses (M ~ 1 Mq for RV, M ~ 0.5 for microlensing) we have placed 
the microlensing point at 8AU, i.e., 3 times the solar-system "snow line" distance. This 
is because microlensing planets are typically detected at 3 times the distance of their own 
systems' snow line (which is of course much closer than 8 AU). 
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Fig. 9.— Planet 



(ICumming et al. 



requencies as a function of mass ratio for mic rolensing (this paper) and RV 



2008 



Mayor et al.ll2009l : I Johnson et al.ll2010l ) detections. As illustrated in 



Fig.lH], the lCumming et al.l ( 120081 ) RV sample and the microlensing sample are consistent with 
each other (despite d i ffereii t freq uencies in thi s figur e) because they are at different distances. 
The ICumming et al.l ( 120081 ) and [Mayor et al.l (120091 ) RV samples are directly comparable on 
this diagram bec ause both are G stars , and they are consistent. However, there is some ten- 
sion between the iJohnson et al.l ( l2010l ) RV measu rement and the microle nsing measurement. 



since both are similar type stars. This is because ICumming et al.l ( l2008l ) and iJohnson et al. 



(120101 ) are "inconsistent" with each other, if one assumes that the frequency of planets as 
a function of mass ratio is inde pendent of host r nass, as we have done in arguing for the 
consistency of microlensing with ICumming et al.l ( 20081). Also sho wn is the (unnormalized) 
slope derived from microlensing observations by ( Sumi et al.lboiol ). 



- 43 - 



-i-> 

> 

S-. 

o 

W 

CD 
C 

a! 



3 



20 



15 - 



10 - 



5 - 







Solar System Analogs 



Neptune 
(0.2) 

Uranus 
(0.3) 



Jupiter 
(5.3) 




1 



Total One Two Three Four 

Planets Planet Planets Planets Planets 

Detected Detected Detected Detected Detected 

(18.3) (5.8) (5.7) (0.3) (<0.1) 



Fig. 10. — Results of Monte Carlo simulation of the anticipated results of our micro lensing 
survey assuming that every star toward the Galactic bulge had a "scaled solar system", as 
specified by Eqs. f lTSj) and f lT6|) . A total of 18 planets would have been detected, including 
6 two-planet systems, compared to the actual detections of 6 planets including 1 two-planet 
system. This seems to indicate that the solar system is overdense in planets, especially 
multiple planets. Also shown are the frequency of various specific combinations. 



